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A source of interchannel interference in certain multichannel FM systems 
is the so-called "frequency 'pulling effect." This effect, which occurs in 
systems using a klystron oscillator, is produced by an impedance mis- 
match between the antenna and the transmission line feeding it. hi this 
paper expressions are developed for the magnitude of the interference wlien 
the speech load is simidated by random noise. 

INTRODUCTION 

In a recent paper' the problem of interchannel interference produced 
by echoes in an FM system was treated. The mathematical development 
in that paper can bo used to caltulatc the diatortion that arises when a 
Klystron oscillator is connected to an antenna through a transmission 
line of appreciable lengtli. 

In the system we study, the composite signal wave (the "baseband 
signal") from a group of carrier telephone channels in frequency division 
multiplex is applied to the repeller of a lOystron and thereby modulates 
the frequency of the Klystron output wave. If the antenna does not 
match the transmission line perfectly, the output frequency is altered 
slightly by an amount proportional to the mismatch. 

This effect, known as "pulhng," results in intermodulation between 
the individual telephone channels. In this study, the composite signal 
will be simulated by a random noise signal of appropriate bandwidth 
and power. It is assumed that some particular message channel is idle; 
i.e., there is no noise energy in the corresponding frequency band (which 
is relatively narrow in comparison with the bandwidth of the composite 
sigiuil). If the system were perfect, no power would be received in this 
idle channel at the output of the FM detector. In the following work, 
the intermodulation noise falling into this chainiel because of the "pull- 
ing effect" will be computed. This leads to "Lewin's integral," so called, 
which is tabulated herein. 
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PULLING EFFECT 

In a perfect FM system the carrier wave can be written 

Eo(t) - A sin [pt + ^(01 (1) 

where A is a constant and the signal is S{t) = d<p/dt = <p'(t), measured 
in radians/second. As mentioned in the Introduction, we assume that 
when the FM oscillator is connected directly to a transmission Hne with a 
shghtly mismatched antenna at the far end, its frequency is changed. 
The reactive component of the input impedance of the line "pulls" the 
frequency of the oscillator to its new value. When the antenna is perfectly 
matched, there is no change in oscillator frequency. 

If the characteristic impedance of the line is Zk and the impedance of 
the antenna is Z^ , the impedance Z looking into the line is 

Zr + Zk tanh P 



Z = Z 



= z 1 + ^^ 



Zk + Zr tanh P 

IP 



1 - pe'-^'- 
where p is the reflection coefficient 

Zh — Zk 



p = 



Z^ + Zk 

and P is the propagation constant of the line. If the loss of the line is 
negligible and the reflection coefficient is small, the input impedance is 
approximately 

Z = Zk[\ + 2p(cos oiT - i sin (j)T)] ohms 

where w is the oscillator frequency in radians per second and T is twice 
the delay of the line. 

It will be observed that the magnitude of the reactive component of Z 
oscillates as the phase angle oiT increases. 

The dependence of the frequency of an oscillator upon the load reac- 
tance has been expressed by earlier workers as a "pulling figure." This 
figure is customarily defined as the difference between the maximum and 
minimum frequencies observed when the load reactance is varied over 
one cycle of its oscillation (the variation being accomplished, say, by 
increasing T). The load is taken to be such that it causes a voltage stand- 
ing wave ratio of 1.5. This corresponds to a reflection coefficient of 0.20 
and 14 db return loss. 

In our work, we assume that the change in frequency is directly pro- 
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portioiial to the reactive component of the input impedance. More pre- 
cisely, \vc assume that the ideal transmitter frequency of p + ^'{t) radi- 
ans/sec is changed by the pulling effect to 

V + ^'(0 + 27rrsiii [Tip + <^'{0)] radians/sec (2) 

where r is given by 

r = 2.5 I p [ X (Pulling Figure in cycles/sec) 

POWER SPECTRUM OF INTERCHANNEL INTERFERENCE 

The distortion produced by the pulling effect is given by the third 
term in (2). This distortion will be denoted by d'{t): 

6'{l) = 27r;- sin {pT + T^'{1)] (3) 

Our problem is to compute the power spectrum of 6'{t). In particular, 
we are interested in the case where the signal ^'(/) represents the compos- 
ite signal wave from a group of carrier telephone channels in frequency 
division multiplex. All of the channels except one are assumed to be busy. 
Although the power spectrum of i^'(0 is zero for frequencies in the idle 
channel, the same is not true for the power spectrum of Q'ii). In fact, 
the interchannel interference (as observed in the idle channel) is given 
hy that portion of the power spectrum of 6'ii) which lies within the idle 
channel. We shall denote the corresponding interchannel interference 
power in the idle channel by tyc(.f)''//' where the idle channel is assumed to 
be of infinitesimal width and to extend from frecjuency / — rf//2 to / + 
dj/2. The function w^ij) will now be computed by using the procedure 
developed in Reference 1. 

The first step is to assume the signal ip'{t) to be a random noise current. 
In order to avoid writing <p' a great many times we shall set tp'it) = S(t), 
where now S(l) stands for the signal. Then the autocorrelation function 
for the distortion 9'(t) is 

Hr{r) = avg[5'(0e'(/+ 7)] 

= {2-Krf avg [sin {Tp + T,p'{t)) sin {Tp + T<p'{t + r))] 

= {2irrf avg [sin {Tp + TS{t)) sin {Tp + TS{t + r))] 

= ^—^ avg [ cos {TS{() - TS(t + r)) (4) 

- cos (2pT + 7'.S'(/) + TS{t + t))] 
- U2irrf (exp [-rVaO) + T'R,(t)] 

- cos (2pT) exp [-T'lUO) - T'Rs(t)]\ 
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where 

Rsir) = [ wsif) cos 2Tfrdf (5) 

•'0 

and uJs(0 is the power spectrum of the applied signal S{t), The last ex- 
pression in (4) follows from the next to the last by analogy with equation 
(1.14) of Reference 1. _ 

The dc component of the distortion e'{t) is its average value e' which 
may be computed from 

f = R,.{^) = ^' [e-^V(l - cos2pr)] (0) 

This follows from (4) since /?s ( « ) = 0. 

The auto-correlation function of the distortion, excluding the dc 
component, is then 

Ji,._j-, = ^ r^V^'lKe'-V- 1) - (e-^'"''''' - 1) oos2pT] (7) 

The interchamiel interference spectrum is 

wl!) = 4: [ Rc{t) cos 27r/T dr (8) 

''0 

where, by analogy with equation (1.22) of Reference 1, 

jlXr) = ^ [e-'^Rme'^'"-''' - T-'Rs(r) - 1) ^^^ 

- (e-^'«-'^' -\-T'Rs{t) -i) cos 2pT] 

As mentioned before, the function wdf) is of interest because 

Pf = w,if) df (10) 

is the average interference power appearing at the receiver in an idle 
channel of width d,f centered on frequency /. 

RATIO OF INTERCHANNEL INTERFERENCE TO SIGNAL POWER 

The average signal power appearing in a busy channel of width df 
centered on the frequency / is 

Ps- ws(!)df (11) 

and hence the ratio of the interchannel interference power to the signal 
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power is 

'i^^An (12) 

p. w,{f) ^ ' 

We now obtain an expression for this ratio on the assumption that the 
random noise signal »S'(0 (which is used to simulate the multichannel 
signal) has the power spectrum 

fn , < / < h 
wsif) = (13) 

lo, / > h 

where Pq is a constant. S(() is measured in radians/sec and Pn/b is meas- 
ured in (radians/sec)^. Pnfb is given by 

Poft ^ ^ = avg y'(t)f = {2iraf 

where a is the rms frequency deviation of the signal measured in cycles/ 
second. According to (5) this signal has the auto-correlation function 

K»(.) = f" Po cos 2./. ,// = P. r?^]' = (2„)' S^ (14) 
Jo L •^"'^ Jo ^ 

whei-e 

U = 2-KfbT 

The interference power spectrum Wcij) corresponding to the Waif) of 
(13) may be obtained by substituting (14) in (9) to get /?c(t) and then 
using (8). The result is 

f e\ A (2xr) _6 \ f, ftu-l sin u 1, -1 • i\ 

WcKj) = 4 — - — e / |(e — hu sm u — 1) 

2 Jo 

(15) 

/ — fi«~l Hill II I I —1 ■ i\ C-, m^ cos (lU , 

— [e -j- hv sm n- — I) cos 2pl\ — ^ - . - du 

2irfb 

where ii is the same as in (14) and we have set 

« = //A h ^ (2iraTf 

This integral may be expressed in terms of Lewin's integral which is 
studied in Appendix III of Reference 1. Thus 

w,(f) = ^'^'^''^f [7(6, a) ~ I( - b, o) cos 2pr](radian/sec)Vcps (16) 
2irfb 

where f{h, a) and I( — b, a) are tabulated for various values of a and b. 
Since we began the problem l)y dealing directly with d'{t) which is a 
radian frequency, rather than 8(1) which is a radian phase, Wc(f) has the 
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Table I — Values of e-''I{b, a) for h > 



b 


e'- 


0.0 


1.000 


0.25 


1.284 


0.5 


1.649 


1.0 


2.718 


2.0 


7.389 


3.0 


20.08 


4.0 


54.60 


5.0 


148.4 


6.0 


403.4 


7.0 


1097. 


8.0 


2981. 



e-^/Cfi, fl) 



a = 



0.000 
0.082 
0.272 
0.761 
1,560 
1.913 
1.974 
1.905 
1.794 
1.680 
1.576 



0.25 



.000 
.072 
,241 
,685 
,440 
,801 
.888 
,844 
.751 
.649 
,552 



0.50 



0.000 
0.062 
0.209 
0.602 
1.291 
1.645 
1.751 
1.731 
1.060 
1,575 
1.492 



0.75 



0.000 
0.052 
0,176 
0,511 
1.117 
1.448 
1.566 
1.571 
1.525 
1.463 
1.398 



1.00 



0.000 
0.042 
0.142 
0.414 
0.919 
1.215 
1.341 
1.372 
1.356 
1.320 
1.277 



Table II — Values of I(b, a) for b < 



1.25 



0.000 
0.031 
0.107 
0.314 
0.713 
0.968 
1.098 
1.153 
1.166 
1.157 
1.138 









Ub 


a) 








a = 


0.25 


0.50 


0.75 


1.0 


1.2S 


0.0 


0.000 


0.000 


0.000 


0.000 


0.000 


0.000 


-0.25 


0.092 


o.oso 


0.068 


0.057 


0.045 


0.034 


-0.5 


0.349 


0.300 


0.254 


0.210 


0.167 


0.125 


-1.0 


1.25 


1.06 


0.885 


0.723 


0.576 


0.432 


-2.0 


4.16 


3.41 


2.76 


2.20 


1.76 


1.34 


-3.0 


8.03 


6.37 


4.97 


3.88 


3.14 


2.46 


-4.0 


12.6 


9.66 


7.23 


5.49 


4.55 


3.74 


-5.0 


17.8 


13.2 


9.40 


6.89 


5.93 


5.19 


-6.0 


23.6 


16.8 


11.4 


8.00 


7.25 


6.85 


-7.0 


30.0 


20.7 


13.1 


8.71 


8.48 


8.78 


-8.0 


37.2 


24.8 


14.5 


8.93 


9.59 


11.0 



dimensions of {radiaiis/se€)Vcps. The signal in the same dimensions is 
Po or (2Traf/fb. Therefore the ratio of the interchannel interference power 
to the signal power is: 



Ps 



2ir 



tl = ±[L] e-\l{h, q) - /( - 6, a) cos2pT] 



(17) 



The quantity e"*/(6, a) for 6 > is tabulated in Table I. The quantity 
/(6, a) for 6 < is given in Table II. These tables, which are also given 
in Reference I, are repeated here for the convenience of the reader. 

When the rms frec^uency deviation a is so small that b = (^2-KaT) 
is small compared to unity, the approximation 



leads to 



I{b, a) ^ ;>V(2 - a)/4 
^ % (2T-m7'-)-(2 - rt)(] - cos 2pr)/2 



(18) 
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When 0- and T are such that 6 » 1 , the approximation 
7(^a)^(6,r/6)^'^exp|^6-^"] 

loads to J; ^{£)' ^ -P [ - ^(2^) ] 

Equation (17), when converted to decibels, breaks down conveniently 
into two terms which may be designated Di and D^ : 

10 log Pi/Ps = D, + D2 

Di = 10 log (r/<jf 

D2 = 10 log ^ e [iQi, a) - /(- b, a) cos 2pT] 
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Fig. 1 — I'lnt of D-_ as a function of aT 



652 THE BELL SYSTEM TECHNICAL JOURNAL, MAY 1957 




-20 



20 



40 60 ao 100 120 

ECHO PHASE, pT, IN DEGREES 



160 



leo 



200 



Fig. 2 — Plot of Z>j as a function of echo phase when f = fb 

The term D2 depends only on the rms deviation a, the round-trip 
echo delay T of the line, the ratio a = f/fb , and the echo phase pT. The 
term Da is plotted in Fig. 1 as a fimction of uT for two channels, one at 
the top and the other at the bottom of the signal band. Since the carrier 
frequency may be expected to be very high, the carrier phase 2pT will 
be a very large number of radians even with very short wave guide runs. 
Hence the curves on Fig. 1 are plotted for the average value of cos 2pT 
which is zero. 

The curves on Fig. 2 show how D2 depends on pT and the parameter 
(tT. In this case the channel is taken at the top of the signal band. When 
pT is an odd multiple of 90 degrees, it turns out that we have even order 
modulation products only; and when pT is a multiple of 180 degrees, 
odd order products only. The curves show that the distortion becomes 
less dependent on the echo phase as the quantity o-T increases. 
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